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g '■ Abstract 

. We explore the possibilities of constructing bulk spacetimes in five dimensions for warped 

braneworld models with a spatially flat Friedmann-Robertson- Walker (FRW) line element on the 
3-brane and with a time-dependent extra dimension. Our first step in this direction involves looking 

. at the status of energy conditions when such a bulk line element is assumed. We check these randi- 



er 



C3 



tions by analysing the relevant inequalities, for specific functional forms (chosen to satisfy certain 
desirable features) of the warp factor, the cosmological scale factor and the extra-dimensional scale 



5h . 

factor. Subsequently, we aim at obtaining solutions with different types of bulk matter sources. 

t^J- . We begin with a general analysis of the solution space of non-singular Randall- Sundrum type bulk 

> 

\& ' models with an exponential warp factor and a chosen equation of state. Thereafter, we focus on 

■ three specific bulk sources - the ordinary scalar field, the Brans-Dicke scalar and the dilaton. In 

. each case, we are able to solve the field equations and obtain desirable solutions for which, we once 

again check the viability of the energy conditions. We also show how one can place branes in the 

o 

bulk using the junction conditions. The issue of resolution of the bulk singularities which appear 



^ ■ in our solutions, using standard methods, is also presented briefly. In summary, we are able to 

demonstrate, that it is indeed possible to construct viable bulk spacetimes for warped cosmologi- 
cal braneworlds with a time-varying extra dimension and with bulk matter satisfying the energy 
conditions. 
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I. INTRODUCTION 



Though yet to be seen in experiments, extra dimensions have been around in theoretical 
constructs for almost a century now. The fact that such a seemingly exotic idea has survived 
for so long seems to make us feel that there is bound to be some amount of advantage in 
having them. The advantage stems from several directions. For instance, one may ask - 
why four dimensions? - a question to which one really does not have an answer. Or, while 
making models of unification (such as superstrings l[) one ends up with a need of more 
dimensions and the subsequent need for compactification. More recently, we have seen how 
the age-old problem -the hierarchy problem-can be solved (obviously not uniquely) with 
extra dimensions 2|. 

The most popular among today's models is the one we mentioned at the end in the last 
paragraph- the so-called warped braneworld or Randall-Sundrum (referred henceforth as 
RS) model. This model assumes a curved spacetime (the bulk) in five dimensions. The four 
dimensional piece however, depends on the extra dimension, a feature unique for this class 
of models (different from the usual Kaluza-Klein 3j). Thus, in such a model we have a 
metric ansatz: 

df ? = e 2f(a) ^2 + ^ + r 2 da 2 

The a = constant (say c7q) section is a four dimensional Minkowski space scaled by the 
factor e 2f ^ =ao) . 

But such a metric ansatz is toy. This is because the metric on the a = constant surface 
must necessarily be cosmological, in order to tally with the universe. Further, there is no 
reason for us to assume r c to be a constant. It could be a function of space and time. Such a 
spacetime dependence effectively makes the scale of the extra dimension, a field, commonly 
known today as the radion. However, here, we restrict ourselves to only a temporal depen- 
dence of r c . Together, these two modifications, i.e. having a cosmological on-brane metric 
and a time dependent extra dimensional scale drive us towards a more realistic braneworld 
model. 

However, the generalisation of the simple RS line element to the more complicated one 
(i.e. with a cosmological scale factor and a time dependent extra dimensional factor) makes 
life difficult for the general relativist to find meaningful exact solutions. Several attempts 
have been made so far. Notable among; them are discussed in . Further, the exis- 
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tence of a time dependent extra dimension in an unwarped scenario has been investigated 
with the aim of finding out whether such a time dependence can drive an accelerating uni- 

il applications including implications for the CMB anisotropies [8j and 
9] have been analysed in recent times. On the other front, warped 
braneworlds with a constant scale for the extra dimension have been studied 
extensively [10| ever since the RS model came into existence. 

Our intentions in the investigations here, is to first figure out whether such solutions 
(one can call them generalised brane-bulk systems) exist and, if they do, how do the matter 
fields, spacetime geometry behave. Our programme is as follows. In the next section, we 
write down the field equations. Then we move on to the issue of energy conditions and try 
to figure out, if it is possible to have viable forms for the unknown metric functions, which 
satisfy them. Then, we look at the possibilities of solutions with matter fields. First, we 
assume a non-singular bulk with an exponential warp factor and discuss the solution space 
for a specific equation of state. Thereafter, we discuss three cases with different types of 
scalar field actions - ordinary and phantom scalar, Brans-Dicke scalar and dilaton. We then 
return to a discussion on the issue of energy conditions for our exact solutions. Finally, we 
discuss the placing of branes using the junction conditions and the question of resolution of 
bulk singularities. We conclude the paper with a summary and some general remarks. 

II. THE LINE ELEMENT ANSATZ AND EINSTEIN'S EQUATIONS 

A general warped line element in five dimensions is given as: 

ds 2 = e 2f{a) g a/3 dx a dx 13 + g u d(J 2 (2.1) 

where g a p can, in principle be any metric and #44 can be a function of space, time and a, 
not necessarily separable. 

The line element we choose to work with is as follows: 



ds 2 = e 2f ^ 



d T ' 2 2j/)2 , 2 -2 



-dr + a 2 it) [ - — + r 2 d8 2 + r 2 sin 



+ rj 2 (t)da 2 (2.2) 



1 — At 2 

In the above, we have chosen a cosmological line element on the a = constant hy- 
persurface. a(t) is the cosmological scale factor and k = 0, ±1 correspond to the usual 
Friedmann-Robertson- Walker (FRW) metric with spacelike sections of zero (R 3 ), positive 



(S 3 ) and negative (H 3 ) curvatures respectively. As defined before, e 2 -^ CT ) is the warp factor. 
The domain of a, which is the extra dimension, could be — oo < a < oo though, as we shall 
see later, we may need to consider finite or semi-infinite domains and/or introduce additional 
branes in order to avoid inevitable bulk singularities which appear in most of our solutions. 
However, unlike the Randall-Sundrum model we now have a time- dependent function r)(t) 
associated with the extra dimension. In general, such a function, when it is dependent on 
all the four coordinates (t, r, 6, <fi), is known as the radion field. It measures the scale of the 
extra dimension at different spatial and temporal locations in the four dimensional world. 
In the two-brane RSI picture, this quantity r\ is known to be the inter-brane distance and 
must have a stable value in order to make sure that the branes do not collapse on to each 
other. Note that in the RS model we always had r](t) as a constant. 

It might be useful at this stage to find out what our expectations are about the nature 
of the unknown functions /(cr), a(t) and i](t). Guided by the RS solution, we expect f(a) 
to be such that the bulk line element is nonsingular w.r.t. the a coordinate, i.e. there are 
no bulk singularities. This, as we will see, is largely impossible unless we postulate it to 
be so. On the other hand, the usual singularity in cosmological time t (or, conformal time) 
should exist. In other words, it is reasonable to assume a big-bang singularity in the four 
dimensional world where a(t) goes to zero at some time t but at all spatial positions r, 9, 
(f) as well as a. Secondly, we would like to have a solution where, following our current 
understanding of cosmology, we have an expanding universe. This will imply a(t) to be 
a monotonically increasing function of t which may have a deceleration or an acceleration 
(negative or positive second derivatives, respectively). On the other hand, we would not 
like to have a large scale for the extra dimension today, which leads to the choice that r](t) 
should be monotonically decreasing in time, but never becoming zero at any finite time. 
This choice is motivated from old ideas in Kaluza-Klein cosmology where, of course extra 
dimensions are always compact. In fact, a growing extra dimension is still fine if it stabilises 
to a certain value at later times. But, as we will see, the solutions found are essentially 
of power law type which, if growing, will keep growing forever. Thus, we choose the extra 
dimension scale to be of the decaying type. We mention that we never use the above choices 
as constraints while solving the Einstein equations. The above choices merely imply that 
we prefer to have an expanding four dimensional universe with a non-compact but shrinking 
(in time) extra dimensional scale. 
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One may argue, following the example in [6j, that an increasing rj(t) is an acceptable 
model essentially because the bulk solution can be transformed to a static, conformally flat 
metric by a suitable coordinate transformation (see [(J). In such a case, a growing extra 
dimensional scale is permissible because it corresponds to the motion of the brane (through 
its embedding) in a static bulk. However, a coordinate transformation of the type used in 

n 

[6|] seems possible only if the Weyl tensor for the bulk geometry vanishes. One can check 
that, for the metric ansatz 12.21 (with k — 0), all the components of the Weyl tensor vanish 
if the following constraint is satisfied, 

£-£ + ^-3 = . (2.3) 
a a z ai] rj 

For solutions of power law type, we may choose, a(t) ~ t Ul and r](t) ~ t U2 . Eq. 12.31 then 
implies 

(v 1 -v 2 )(v 2 -l) = 0. (2.4) 

v\ = V2 means a trivial solution which results in the same functional form for both the scale 
factors and v-i — 1 is essentially what was found in {(|. It may be noted that, for the scale 
factors obtained as solutions in our examples, the above constraint relation does not hold 
and thus, the bulk geometries are not conformally flat. 

Furthermore, following the discussion in 0], a moving brane arises when we make a gauge 
choice in which the metric on the a — t section is manifestly conformally flat. Fixing the 
brane implies that we do not retain the full gauge freedom. Thus, the growing/decaying 
nature of the extra dimension is linked with the motion of the brane in the bulk. By making 
a choice of decaying extra dimensions, we therefore, do not eliminate the possibility of viable 
models with growing extra dimensions though our preference (and the obtained solutions) 
is based on some logic- essentially linked with old Kaluza-Klein ideas and the issue of a 
stabilised extra dimensional scale at later times. We mention that such solutions with a 
decaying extra dimensional scale have been discussed earlier, for example in 4| itself. 

Thus, our central question now is: do the Einstein field equations with a matter source 
have such kind of solution that we seek? To answer this query, we must write down these 
equations. We first choose to write down the Einstein tensor and perform our analysis by 



assuming it to be the required stress energy for specific choices of the unknown 
This, in the older literature of general relativity is known as Synge's g-method 



unctions. 



ll|. The 



other way, is to actually solve the Einstein's equations with some specific choices of matter 



stress energy (eg. different types of scalar fields residing in the bulk) - this, following Synge 
is the T-method. 

The nonzero components of Einstein tensor (in the frame basis) for the metric 12.21 with 
k = (i.e. a spatially flat cosmological brane) are, 

Goo = e- 2S U~ 2 + S-t) - ~ 2 (3f + 6/' 2 ) (2.5) 
\ a z a rj J 7] z v ' 

G aa = -e-' { 2* + % + 2 a -1 + 5) + -i (3/» + 6/' 2 ) (2.6) 
\ a a z ar\ r\ J rj z v 7 

2f / a a 2 \ 6f' 2 , . 
G 44 = -e~ 2/ 3- + 3^ + -V 2.7 
y a cry rj z 

and G 04 = ^f'e~ f (2.8) 

where, in the above, G aa refers Gn, G 2 2, G33, the three spatial indices in the 4D FRW metric. 
A dot denotes differentiation w.r.t time while a prime indicates differentiation w.r.t. a. 
Defining new variables, ^ = x, 2 = y, we have, 

G 00 = e"^(3x 2 + 3x2/) - ~ 2 {?>f" + 6f' 2 ) 

T] z 

G aa = -e^{2x + 3x 2 + 2xy + y + y 2 ) + 1(3/" + 6f 2 ) 

rj z 

C 44 = - e ~ 2 f (3i + 6x 2 ) + 4t6/ /2 

and G 04 = 3t//' (2.9) 

77 

Some interesting features in the structure of these equations may be noted here. 

• The presence of the flux term emerging out of G04 which is zero only when y = 2 is 
zero and/or when the geometry is unwarped (i.e. / is a constant). 

• The separability of functions of time and a which enables us to make an attempt 
towards constructing solutions by equating coefficients of e~ 2 f and -5 on both sides of 
the Einstein's equations with some matter source. 

• A vacuum solution (Gjj = 0) turns out to be either un-warped or with a time indepen- 
dent extra dimension. Also, with a bulk cosmological constant, (i.e. Gij + Agij = 0), 
a non-trivial solution cannot be found. 
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III. THE STATUS OF ENERGY CONDITIONS 



Given the above expressions for the Einstein tensors, we can now ask the question: is it 
possible to have geometries, i.e. specific functional forms ofa(t), r](t) and f(a), with matter 
satisfying the weak or null energy conditions? To do this analysis we need to know about the 
energy conditions in a bit more detail, largely because certain nontrivialities arise because 
of a nonzero G04 term (equivalently, a nonzero T 4 must also be there) being present. 



A. Energy conditions: weak and null 



The Weak Energy Condition is given by, 

GuU 1 ^ > (3.1) 

where U 1 is a nonspacelike vector and a similar inequality defines then Null Energy Condition 
if U 1 is a null vector. 

One can write down the energy conditions in terms of the eigenvalues of the energy- 
momentum tensor (which must be real). If Ao denotes the eigenvalue corresponding to the 
timelike eigenvector, the weak energy condition is equivalent to the following simple relations 
amongst the eigenvalues [121 ] . 

-A >0 and - A + A Q > where (a = 1,2,3,4). (3.2) 

In effect, we find the inequalities by diagonalising the energy-momentum tensor. The eigen- 
values of the energy-momentum tensor are the roots of the equation 



\Gu — Xgu\ 







(3.3) 



which reduces to, 



where, 



p+ A 
Pl - A 


Q 





Pi- X 
pi — A 




3/" + 6f 2 
rf 



a 



(1 




P2 - A 
arj 



+ e~ 2f 3— + 3- 

V or an 



(3.4) 



(3.5) 
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Pl = JL ( 3/ » + er) - e-V (2° + 4 + 2*5 + 5) , (3.6) 
v y \ a or at] rj J 



6f' 2 2f ( a a 2 \ 
^ = -^"^ 2/ 3- + 3- (3.7) 



and q = e~ 2f . (3.8) 



if 

Here, and henceforth, we shall absorb the factor k in the Einstein equation Gjj = kTjj 
through appropriate scaling redefinitions. The five eigenvalues from Eq. 13.41 assume the 
following forms, 

Ao = \{~P + P2 - v /(p + p 2 ) 2 -4g2) (3.9) 

Ai = A 2 = A 3 = pi (3.10) 

1 



A 4 = -(-p + p 2 + yj( p + p 2 f - 4g 2 ) (3.11) 
The fact that the above eigenvalues have to be real, leads to the following requirment 

(p + p 2 ) 2 -4g 2 >0. (3.12) 

This may restrict the allowed domains of t and a. As the eigenvector corresponding for the 
eigenvalue An is timelike or null, the effective inequalities for Weak Energy Condition turns 



out to be 



131], 



Ft = p-p 2 + V(P + P2) 2 - V > 0, (3.13) 



F 2 = p + 2pi - p 2 + J(p + P2) 2 ~ M 2 > 0, (3.14) 



and F 3 = v / (p + p 2 ) 2 ~V > 0. (3.15) 
FigfS] shows status of the above three inequalities for two models in the two columns, 

where f(a) = —b\og(cosh a) represents a thick brane, 

a(t) ~ e Ht a de-Sitter brane (3.16) 
and r](t) = c + ee~ m 

(where t and a run from to 00 and —00 to 00 respectively) with e = 1 and e = — 1 
respectively (variation of rj(t) with two different values of e is shown in Fig{T]). Figj3] shows 
the variations of the same functions (i.e. ^1,2,3) for 

a(t) ~ \ft a radiative brane 
and n(t) = c + ee-P H{t ~ to) (3.17) 

8 




FIG. 2: Status of the inequalities for a de-Sitter brane for following set of parameter values: b = |, 
c = 2, (3 = 3, H = 1, e = 1 (left column) and e = — 1 (right column). 

where t runs from t (= 1 say) to oo and the warp factor is the same as in Eq. 13.161 

We find that for the above two models, all the inequalities are satisfied at least with the 
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FIG. 3: Status of the inequalities for a radiative brane for following set of parameter values: b = 2, 
c = 6, (5 = 2.4, H = 1, to = 1, e = 1 (left column) and e = — 1 (right column). 

specific choices of parameters we have made. It has been checked that these inequalities can 
be satisfied in the presence of a decaying warp factor. However, for a growing warp factor 
(i.e. for a negative b) the inequalities are violated at least in some spacetime region. 

IV. EXACT ANALYTIC SOLUTIONS IN THE BULK 

In this section, we discuss possible exact solutions in the bulk. We first assume an 
exponential warp factor and discuss some specific consequences for a given equation of 
state. In subsequent sub-sections, we look at exact solutions with different types of scalar 
fields - an ordinary massless scalar, the Brans-Dicke scalar and finally a dilatonic scalar. 

A. Models with an exponential (Randall Sundrum type) warp factor 

In the Randall-Sundrum scenario, the warp factor is of the form f(a) = — b\a\. Using the 
functional form f(a) = —ba (or bo) and restricting ourselves to the a > (or a < 0) region 
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of the extra dimension, we now investigate some special cases by looking at the Einstein 
tensors and the resulting required matter stress-energy. To make things more quantitative, 
let us first write down the Einstein tensors with the above-mentioned warp factor. 

6b 2 

Goo = e 2b °(3x 2 + 3xy)-^, 

6b 2 

G aa = -e 2ba {2x + 3x 2 + 2xy + y + y 2 ) + — , 

6b 2 

G 44 = _ e 2fcCT (3x + 6x 2 ) + ^- r 

i] z 

e ba 

and G*o4 = 3yb — (4.1) 

V 

(a) Let us first look at the situation where we have 

Goo — —G aa = — G44. (4-2) 

Since the above relation holds for the factors associated with the -4 terms in the Einstein 
tensors, conditions on x, y and its derivatives emerge when we impose identical conditions 
for the factors associated with the e 2bcr terms. These conditions, after some elementary 
algebra lead to: 

x = —x 2 + xy, 

y = 2x 2 — xy — y 2 . (4.3) 

Eq. 14.31 is an autonomous dynamical system of nonlinear first order differential equations. 
The general feature of all the solutions can be determined through a solution space analysis 



14j of the above system. In this case, the line x = y is a critical curve, i.e. every point on this 
line (except the origin) is a non- isolated fixed point. The eigenvalues of the Jacobian matrix 
(of the linearised system) at any such fixed point (x*, y*) are (0, — 4a:*), which suggests, for 
x* > fixed points are neutrally stable whereas for x* < they are neutrally unstable. It 
is difficult to make any conclusive comments, analytically, on the behaviour of phase space 
trajectories near those points, because such borderline cases are very sensitive to nonlinear 
terms in the equations. But the phase portrait (FigHJ) indeed confirms our conclusions about 
the nature of the fixed points. Note that, the (0, 0) point in the phase space is irrelevant 
as the Jacobian itself becomes singular and therefore it does not fall in any category of 
fixed points. From a physical point of view, the origin represents a static universe. We 
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will see that no trajectory reaches that state in finite time. This provides a justification for 
the origin not representing a viable solution. Now, we turn our attention to the physically 
most preferable region of the solution space. This is the below-right quadrant, where x is 
positive and y is negative. There is one such curve, x = —y, which is distinct from every 
other. It approaches the origin (0,0), i.e. towards a static universe, as t — > oo. Any (x,y) 
combination that is located above that line flows toward the line x = y, resulting in the 
same scale factor for both kind of spatial dimensions, at t — oo. On the other hand, flows 
located below the curve x = —y, tend toward a static brane (x = 0) with an extra dimension 
whose size shrinks, with ever increasing rate (y — > — oo), towards zero as t — > oo. 
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FIG. 4: The solution space for the dynamical system in Eq. 4.3. The solid lines indicate different 
flows in the solution space of {x,y} for the following different initial conditions: {x(0),y(0)} = 
(0.5, -0.5), (0.47, -0.55), (-0.05, 0.8) and (-0.11, 0.11), for t running from to 3.5. 

Following the phase portrait, we assume y = nx, which implies, for consistency, n 2 = 1 
or n = ±1. The n=l case is trivial because it leads to constant a(t) and hence, constant 
rj(t). With n = — 1 we obtain, a(t) ~ and r](t) ~ t~^. This specific solution represents a 
radiative brane while the extra dimesion is indeed decaying (in accordance with our prefer- 
ence). For these forms of the a(t) and r)(t) the stress energy is remarkably simple because, 
the coefficients of e 2ba in the Einstein tensors are all zero. Therefore, we just have: 

Qb 2 3b 
Too — — — — —T aa — T 44 ,T 4 — — . (4.4) 
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(b) We now move on to another case where we choose a(t) = e 131 and n(t) = e~ yt with 
/?, 7 > while f(a) remains the same (i.e. f(a) oc a). The Einstein tensors now take the 
form 

e -2 7 t ' 



Goo = e 2bCT (3/? 2 - 3/? 7 ) - 



G aa = -e 2feCT (3/3 2 -2/3 7 + 7 2 ) + 
G 44 = -e 2b °(Q(3 2 ) + 



6b 



.2 



2 7 i ' 

6b 2 

,-2 7 t ' 



6(7 

G 04 = -37^- (4.5) 

Notice that if /9 = —7, we end up with the same relations between the Einstein tensors as 
in (a) above, though this involves growing factors for both the cosmological scale and the 
extra dimension scale. In addition, we also note that with 7 = one obtains a de-Sitter (or 
anti de-Sitter) brane with a constant scale for the extra dimension. 

We now move on towards obtaining solutions in the true sense by assuming specific forms 
of the bulk energy momentum tensor. 

B. Bulk ordinary scalar 

The simplest choice for a bulk energy momentum tensor is that of an ordinary, massless 
scalar field, for which we have: 

Tjf lar = d^dj<\> - l - gi jd K (t>d K <P (4.6) 

with components, such as 

Too = ^-0 2 + — 2 = T 44 , 

T aa = —<P 2 -^ and T 04 = #'— . (4.7) 

We now need to equate the above with the Einstein tensors and obtain solutions. To keep 
things simple, let us assume 

<K*,<r) = <Mt)+<M<7). (4-8) 

From the equation for the 04 component, we note that <p\ = my and <p' 2 = ^/'. So, the form 
of (j)(t, a) will be 

3 

(j>(t, a) = m log[77(*)] + — f(a) + const. (4.9) 
13 



Using the relations between the coefficients of \ on both sides gives: 

f(a) = i log(<r - (Tq) + constant ; f' 2 = 12(f)' 2 ' \ (4.10) 
In the same way, using the relations between the coefficients of e -2 ^ one gets: 

y = -y 2 - 3xy. (4.11) 

We now further assume y = nx which is consistent, for all n, with the above two equations 
for x and y. Therefore, solving the above two equations one finds 

1 n 

a(t) ~ (*-ti)*rs ; rj(t) ~ (t-t 2 )^, (4.12) 

where t\ and t% are integration constants. However, overall consistency requires that m 2 = | 
and n = 4±2v^6- For n = 4 + 2a/6, both a(t) and rj(t) are growing functions of time, which is 
a feature not desirable. On the other hand, for n = 4 — 2y/E, we have the proper behaviour 
for a(t) and n(t). The values of the exponents for a(t) and n(t) are quoted in Table 1. 
It is worth mentioning here that solutions of Eq. 14.111 are constrained by a consistency 
requirement of type x oc y, with a specific proportionality constant, which leaves us with 
solutions very few in number with respect to what we obtained in the case discussed earlier. 
Thus, a solution space analysis of Eq. 14.111 becomes irrelevant in this case. In fact, similar 
conclusions apply to the other two cases analysed below. 

The case of the phantom scalar (with a negative kinetic energy) has also been investigated. 
Using the same methods, we find that solutions do not exist, primarily because, we end up 
having 0' 2 = — 12/3 2 /' 2 , which is impossible, unless both sides are identically zero. 



C. Bulk Brans Dicke scalar 

Brans-Dicke theory [16] is well-known as an alternative theory of gravity where a scalar 
field is assumed to be responsible for generating the gravitational constant G. Though 
experimentally almost ruled out, it serves as a useful model and, as mentioned later (see next 
section on dilaton gravity), it has also reappeared in various contexts in recent times. There 
have also been a few attempts [l?j, [3] towards constructing warped braneworld models in 
five-dimensional Brans-Dicke theory. 
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The action for Brans-Dicke gravity (in five dimensions and without a potential) where 
the Brans-Dicke scalar and the metric are the basic fields, is given as, 



S 



f ofxJ—q 

16tt J v y 



(j)R 



lK 



matter 



(4.13) 



In the above, we have assumed additional matter other than the scalar field itself (which, 
in a sense is not really matter, as such). The resulting Einstein's equations and the scalar 
field equation are given as, 



G 



1.1 



-T L U + To 



<P,i4>,j - ^9ij4>,k(P k 



+ 



;J;J ~ 9lJ 





□ ( 



□0 (V0) 2 

Contracting Eq. 14. 141 and substituting in Eq. 14.151 we get 

8tt 



-T 



4 + 3cj 

where T denotes the trace of the matter energy momentum tensor. 



(4.14) 
(4.15) 

(4.16) 



1. Solutions with perfect fluid bulk matter 

Let us choose the bulk matter energy momentum tensor to be that of a perfect fluid, 
Tjj = diag(p, p a , p a , p a , Pa) , with a vanishing trace 



T = -p + 3p a + Pa = 0. 



(4.17) 



Thus the scalar field equation essentially becomes 



□ c 



-2/ 



+ s5 + 2U 

V a Tj J 



4 W + */' 



Then assuming 4>{t,a) = 0i(£)0 2 (cr), the Einstein's equations yield, 



G 



00 



-2/ / 01 , U< Pt 



i,'2 



0! ' 2 0f/ r/ 2 V 2 ^ 2 



(4.18) 
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e-f 

and G04 = 

V 



0102 01 02 



(4.19) 

Using the expressions for the Einstein tensor components, as given in Eq. 12.91 the off- 
diagonal term in the Einstein equations leads to the following possibility, 

— = my and — = nf where m = ? U . (4.20) 

01 02 f-1 ^ J 

Thus, we have, 

0(t,a) ~77 m (t)e" /(<T) . (4.21) 

Now, dividing Eq. 14.181 by 0, separately equating the coefficients of e~ 2 ^ and \ to zero 
and using fl~20l we get, 

y + 3xy + (m + l)y 2 = (4.22) 
and f" + {n + A)f' 2 = (4.23) 

Similarly, using diagonal components of Einstein equations and Eq. 14.22^ tracelessness 
of T u (Eq. KT7\f leads to, 



x + 2x- Ay 2 = 0, where —J (4.24) 

-^-n 2 -n 2 -4n + 30J f' 2 = (4.25) 

Consistency requirement between Eq. 14.231 and Eq. 14.251 gives rise to the following condition, 



4±2v/4T3^ 3 tA 
n± = which also implies m± = . (4.26) 

Then Eq. 14.231 results in 

f(a) = const. + — — \og(a - a ) . (4.27) 
4 + n 

It may be noted that, the factor is always positive for n = n_, but it can be negative as 
well for n = n + , which implies both growing and decaying warp factor solutions are possible. 

It is evident that Eq. 14.221 and Eq. 14.241 constitute a dynamical system of their own for 
every possible value of u. At the critical point (0,0), the Jacobian of the linearised system 
vanishes. Therefore, let us directly look at the phase portraits (Fig. [5]) of this system for 
m = m + (left plot) and m = m_ (right plot) at, for example, u — 1. Cosmologically 
viable features are apparent in the x-positive (expanding brane) and y-negative (decaying 
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FIG. 5: x — y solution spaces at to = 1 with m = m + (left plot) and m = m_ (right plot). Four 
trajectories in both the figures are drawn for same initial conditions, which are - {x(0) = 0.5, y(0) = 
0.5}, {x(0) = 0.5, y(0) = -0.5}, {x(0) = -0.3, y(0) = 0.2} and {x(0) = -0.2, y(0) = -0.2}, while 
t runs from to 1. 

extra dimension) quadrant of both the phase spaces. In the left-side plot, one class of 
trajectories flow towards the centre, as time increases. Another class of trajectories flow 
toward y = — oo (extra dimension decaying at an ever increasing rate) solution as t — > oo, 
while the x-coordinate increases very less. Now, x becoming almost constant, in turn implies 
an almost exponential growth of a{t) or brane inflation. In this case, this can mean that the 
brane expands very rapidly and forever - a big rip. 

On the other hand, in the right-side plot, one class of x-positive and y- negative solutions 
flow towards the centre , which represents a static universe and another class of trajectories 
flow toward x = — oo and y = — oo (whole of the space is singular at this point) while 
crossing the y-axis (x = 0) very slowly. This means the expansion of the brane first slows 
down to staticity and then starts contracting - a big crunch. Another crossover takes place 
from x- negative and y-positive quadrant to x-positive and y-positive quadrant, for certain 
initial conditions, and then these trajectories flow towards the centre. Some part of these 
trajectories become parallel to y-axis, i.e. x become constant, for a very short span of time, 
which can cause brane inflation, though the extra dimension is also of the growing type. 
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2. The vacuum solution 



Now, let us assume 



y = hx, 



(4.28) 



i.e. we are considering the straight lines in the phase spaces at all values of u. Then the 
consistency requirement between Eq. 14.221 and Eq. 14.241 leads to the following condition 



Ah 2 + (m + l)h 2 + 1 = 



m + 1) ± yJ(m+l) 2 -AA 
2A ' 



(4.29) 



When the abovementioned constraint is satisfied, our job reduces to solving only one inde- 
pendent equation 

x = (Ah 2 - 2)x. (4.30) 
Integrating further, we eventualliy get 



a(t) ~ f*- Ah ' 2 



and 



i](t) ~ t 2 ~ Ah ' A . 



(4.31) 



The most interesting feature of this specific set of solutions is that, these are actually vacuum 
solutions, i.e. one can easily check that all the matter energy momentum tensor components 
identically vanish when Eq. 14.291 is satisfied. The following figure gives a view of the 
exponents in a(t) (yi = 2 __^ fe a ) and rj(t) (u 2 = 2 _ ) as functions of u. It can be clearly 



h = h + 



h = h 




0.4 
0.2 



-0.2 
-0 . 4 



Vl 

y 2 



10 



FIG. 6: Variation of the exponents, V\ and V2, of the scale factors a(t) and n(t) with respect to to 
for h = h + and h = h_ are shown in the left and right-side plots respectively. 

seen that all possible combinations of v\ and u 2 (+ve+ve, +ve-ve, -ve+ve and -ve-ve) are 
there at different values of uo. In fact, the variations seem to suggest that an expanding 
brane mostly suits a decaying extra dimension. 
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3. The radiative solution 



One specific case, where an analytic solution is possible, is when A = 0, i.e. for to = 1.586, 
we have 



x 



-2x 2 



-3xy - 2.26y 2 . 



(4.32) 
(4.33) 



Which implies, 



1 



and 



x 



y 



2t-c x 



1 



(4.34) 
(4.35) 



(2t - c x ){c 2 \/2t - c x - 2.26) ' 

where c\ and c 2 are integration constants, whose different values will span the whole solution 
space of x and y for the abovementioned value of u. It may also be noted that, with c 2 = 0, 
we have the vacuum solution for uj = 1.586 (in which case the extra dimension is essentially 
decaying again). Further, integrating Eq. 14.341 and Eq. 14.351 we get, 



and 



a{t) ~ (2t-d)5 
. (2.26 - c 2 V2t^) 0Ai 

7]{t) 



(4.36) 
(4.37) 



(2t - Cl ) - 22 

The above solution, in fact, represents a radiative brane while the nature of the extra di- 
mension depends on c\ and c 2 - Fig. [7| shows how the nonzero components of matter energy 
momentum tensor, p and p a (p a vanishes in this case), behave as function t and a with 
ci = and c 2 = 1. 





FIG. 7: Behaviour of the nonzero components of energy momentum tensor for bulk matter with 
radiative brane (using m = m + and n = n + ) with c± = and c 2 = 1. 
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4- Solutions in Brans-Dicke frame using Einstein frame scalar field solutions 



It is well known that the action of Brans-Dicke theory, as stated in the previous section, 
can be converted into that of canonical Einstein gravity coupled to a massless scalar field. 
We recall below how this is done. First, let us rewrite the action for Brans Dicke theory 
replacing <fi by e~ 2 ^ and gu by gu . Subsequently, defining a conformally related metric 

~g I j = tt 2 {a,t)g IJ (4.38) 

and choosing f2 = we find that the Brans-Dicke action (with no additional matter fields) 
goes over to: 



S = J Sx^~g 



(4 + 3w)0 JJ Vj0Vj0 



(4.39) 



Further, defining: 



$ = -v/4 + 3u;0 (4.40) 
o 



we can convert the action above into that of Einstein gravity coupled to an ordinary scalar 
field. 

This equivalence can now be used to construct new solutions in Brans-Dicke theory by 
making use of the solutions with an ordinary scalar field discussed earlier in Section IV 
B. The main point here is that the solutions with an ordinary scalar are also solutions of 
Brans-Dicke theory in the Einstein canonical frame. How do these solutions look like in the 
Brans-Dicke frame? We look at this aspect now. 

Following the earlier ordinary scalar field analysis we choose: 

3 

& = mlnri-\ / (4.41) 

m 

with a(t) ~ (t-to)^, rj(t) ~ (t-t )™+* and f{a) = \ \n(a-a ). We also had m 2 = f and 
n = 4 ± 2y/6. 

The metric in the Brans-Dicke frame is related to that in the Einstein frame by an overall 
conformal factor f2 2 = e^. Using the relation between and $ we find that the conformally 
related line element becomes: 

ds 2 = e 2Mai) [-dr 2 + al(r)\dx\ 2 ] + vK^daf (4.42) 

where: 

r r ( * ~ t0)P+1 " « « ^ ~ ^ (AA^ 

T - T0= p+1 ' ai - ai= q+l (443) 
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p 



ran 



4mJ3(4 + 3u;) 



(4.44) 



and 



/i(<ri) = const. + Aq + 1 log(a- x - q-i) , 
4g + 4 

p(n + 3) + l 

ai(r) ~ (r - r ) (p+ 1 )("+ 3 ) , 

p(n + 3) + n 

r/i(r) ~ (r - To ) (p +1 >< n+3 > . 



(4.45) 

(4.46) 
(4.47) 



We now illustrate the nature of the coefficient in fx and the exponents in ci\ and n\ appearing 
in the above solutions as functions of u for the four possible combinations of (m,n) through 
Fig. [BJ It is interesting to note that for set(i) and set(iii) growing ai(r) - decaying rji(r) 



(i) 




I 0.5 



-1 1 """2 3 A 5- 

(iii) 



~2"' 3 4 5 



-1 

-0.25 
-0 . 5 
-fl". 75 
' -1 



(ii) 



12 3 



(iv) 



1 2 3 4 5 



FIG. 8: Variations of the coefficient 



4g+l 

4y+'l 



p(n+3)+l 
(p+l)(n+3) 



in the extra 



in /i(o"i) (continuous line), the exponent 
in the cosmological scale factor a\{r) (dotted lines) and the exponent 
dimensional scale factor tji(t) (dashed lines) w.r.t. oj are shown with following different (m, n) 
combinations; set(i): h§, 4 + 2y/&), set(ii): 4 - 2^) , set(iii): (=^, 4 + 2^) and set(iv): 



V3 



2V6 • 



combination does not exist whereas for set (ii) and set(iv) both decaying and growing warp 
factor solutions exist with desired evolutions for ai(r) and ^i(t). 
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D. Bulk dilaton scalar 

Low energy effective string theory gives rise to Einstein-like equations through the con- 
ditions that the /^-functions of the string c-model are equal to zero The Einstein like 
equations have additional terms involving the dilaton (a scalar), third rank antisymmetric 
tensor and other fields (Maxwell and moduli fields) which arise out of the method of com- 
pactification. Dilaton gravity involves only the dilaton and the metric field. It is different 
from ordinary scalar field theory and is equivalent to Brans-Dicke theory under a special 
choice of the u parameter (being set equal to -1). The action for dilaton gravity [15( is given 
as: 

S = J d 5 x^ge- 2 ^(R + Ad K (f)d K (p). (4.48) 
This gives the following field equations, 

Gjj = -2 V/V A + guV 2 <P, (4.49) 
</"V/Vj0 = 2g IJ V I <f>V J <f>. (4.50) 



1. A solution in the string frame 

The terms in the R. H. S. of the above equation can be clubbed together to yield an 
effective energy momentum tensor. The effective matter stress energy is therefore given as: 

ie. Too = e~ 2/ (-0 + (3x + y)j>) - \(<//' + 2f<j>'), 

T aa = -e" 2/ (0+(x + ?/ )0) + ^((/) // + 2/V / ), 

T 44 = - e - 2 /(0+(3x-y)0) + ^(-<f + 4/V) 

r] z 

and T 04 = { r 2d t d a <t> + 2f'<l> + 2y<tJ). (4.51) 

V 

Let us assume, as before, 

<f>(t,a) =Mt)+M*)- ( 4 - 52 ) 
Then the off-diagonal term in the Einstein equations gives us the following constraints, 

3 — 2n 

0i = my and <f> 2 = nf where m = . (4.53) 
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So, essentially we have, 



(f)(t, a)—m \og[r)(t)} + nf(a) + const. (4.54) 
Now, equating the coefficients of in the both sides of the Einstein's equations we get, 

/» + ^V'V) = (4.55) 
and (m,n)=^l) or (~,3). (4.56) 



On the other hand, equating the coefficient of \ in the scalar field equation gives: 

+ = (4.57) 



Thus the only allowed value for (m,n) is m — \ and n — 1 which gives: 

f(a) — const.-\ — log((T — (To), a growing warp factor, (4.58) 

Finally, equating the the coefficients of e~ 2f in the Einstein equations (with m — \), we 

get, 

2 3 1 1 2 
3x + -xy + -y - -y =0, 

2 3 1 1 2 
2x + 3x + -xy + -y+ -y =0 

13 1 

and 3x + Qx 2 - -y - -xy + -y 2 = 0. (4.59) 
The three coupled equations above give rise to the following algebraic constraint, 

3a; 2 _ - y 2 = o. (4.60) 
2 



which yields 



y = ±V6x, - + 2- = and 5 + ( ±J* _ 1 ) t = , (4.61) 



with solutions as, 

a(t) ~ (t - to) ^ and r/(i) ~ (t - t ) ±v/ ^- (4.62) 

Note that the coefficient of e~ 2f in the scalar field equation gives y + 3xy = 0, which is 
automatically satisfied by the abovementioned solution. We discard the positive exponent 
solution for rj, on physical grounds discussed before. 
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2. Solutions in string frame using Einstein- scalar solutions in Einstein frame 



It has been noted earlier that for uj = — 1, Brans-Dicke theory gives rise to dilaton gravity. 
Using the solutions for Einstein gravity coupled to an ordinary scalar we now construct string 
(Brans-Dicke) frame solutions in dilaton gravity, following the discussion presented in IV 
B. 



Substituting u — — 1 in Eqns. I4.45H4.47I we get the following forms for ai(r), 771 (r) and 
AM- 



/i(cri) = const. + q + ] log(<7i - 0-1) 
4g + 4 

p(ra+3) + l 

0i (t) = [(p + l)(r - t ] (p+i)(«+s) 

p(n + 3) + n 

771 (t) = [(p + l)(r - r )] cp+dc+3) 



(4.63) 

(4.64) 
(4.65) 



with p 



and m = n = 4 ± 2^6 



%/3(n+3)' ^ " I'" '"" L "' " ~ :? ' 



The four sets for ai, 771 and fi (correspomding to four different combination of (m,n)) 
yield two acceptable solutions - one with a decaying warp factor and another with a growing 
warp factor (this is in fact the same solution as given by Eq. I4.62j) . Figj9] shows the nature 
of these solutions in presence of bulk dilaton field. 



fi(o-i) 




3.5 1 1.5 2 2.5 3 



f l(o-l) 



ai(r) 

: : ! 



r\ 2 3 



0.5 1 1.5 2 2.5 3 



0.5 1 1.5 2 2.5 3 



FIG. 9: Nature of the possible combinations of ai(r), ryi(r) and /i(o"i) found in presence of a bulk 
dilaton field. In the top row: (m,n) = (-^,4 ± 2\/6) and in the bottom row: (m,n) = (^-^,4 ± 



2\/6). The exponents in ai(r) are same for 



m,nj 



;^,4 + 2\/6) and 



m,n) 



4 -2y/6). 



The explicit solutions (those which are physically meaningful) are displayed in Table I. 
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V. STATUS OF ENERGY CONDITIONS FOR THE ABOVE SOLUTIONS 



We now focus our attention on analysing the energy conditions some of the above solu- 
tions. 



Bulk Field 


/(<t) = const. + 


a(t) ~ 


v(t) ~ 


Fx > 


F 2 >0 


F 3 > 


Ordinary Scalar 


0.251og(cr-CT ) 


it - to) ' 48 


(i - *o)- a43 




X 


v 7 


Brans-Dicke (for to — 1) 


0.371og(cr — o"o) 


it - to) ' 43 


(i - *o)-°' 55 


V 


X 


v 7 


0.075 log(o--CTo) 


it - 1 ) ' 515 


(f - to)- ' 32 


V 


V 


v 7 


Dilaton Scalar (uj — —1) 


0.51og(<7 - cr ) 


(t-to) - 33 


(t-to)" - 82 


X 


X 


v 7 


-0.51og((T - (7 ) 


(t - t f 57 


(t - to)" - 18 


v 7 


v 7 


v 7 



TABLE I: The above table shows few of the typical solutions, found in Section IV, mentioning 



their desirability as viable models and the status of the energy condition inequalities (Eqs. 13.131 

EES EE]). 

In Table UJ we have listed the 'good' solutions (only those which have a combination of 
growing cosmological scale and a decaying extra dimensional scale) found in the previous 
sections. For the Brans-Dicke case ,we have two solutions for u — 1. The solutions with 
a dilaton field are essentially derived from general solutions for the Brans-Dicke case by 
just equating uj — — 1. There are solutions with both decaying and growing warp factors. 
Notably, all solutions have singularities at a finite value of t or a. To investigate the status 
of WEC for these solutions, we have plotted the inequality functions. It is found that 
the 3rd and 5th set of solutions satisfy all the inequalities (Fig iTOl . Here, the statement 
inequality is satisfied means that it is satisfied in an entire semi-infinite spacetime region 
from {(t — t c ,o~ — a c ) to (t = oo,a = oo)}, where t c and a c are suitably chosen lower 
bounds (these bounds essentially comes from the constraint inequality relation 13. L2[) on the 
allowed domain for t and a (energy condition satisfying solutions have to be defined in these 
domains only). The statement not satisfied means that no semi-infinite spacetime region 
can be found where energy conditions are satisfied. 
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FIG. 10: Status of the inequalities for the two "good" solutions as given in the previous table for 
to = and ctq = with suitably chosen t c and a c . 

VI. PLACING THE BRANE AND SINGULARITY RESOLUTION 

Till now we have essentially found bulk solutions. To have a proper brane-bulk system, 
we need to place a four dimensional hypersurface in the bulk. One can do this by replacing 
the argument of the logarithm (appearing in every solution of /(cr)), i.e. a, by \o~\. The 
contribution of this modulus function will appear as a Dirac delta function peaked at some 
cr& (which is the location of the brane), in the left hand side of Einstein equations. To justify 
this, the energy momentum tensor on the right hand side of Einstein equations must carry 
a delta-peaked term too. This can be done by adding an extra term multiplied with a delta 
function in the bulk Lagrangian, which also peaks at the location where we intend to place 
the thin brane. 



A. Junction conditions 



Let us take the contribution of the brane to the Einstein-Hilbert action to be, 

Sb = - [ d*xy/=& A(0) = - / d 5 x^9~ 5 X(<j>) 6{a Zf b) . (6.1) 
J J Jg^ 



Here g§ denotes the determinant of the full five-dimensional bulk metric and g^ is the deter- 
minant of the induced metric on the brane. Now, we can determine the jump in f'(cr) (or 
4>'{o-)) across cases of the ordinary 

scalar field or the Brans-Dicke field coupled to gravity. 
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1. Junction conditions with ordinary scalar field 



Integrating the scalar field equation at the location of the brane, we get, 



Then from Einstein equations we have, 



T) dX 



If] - -2a 



<r=a b 



a=a b 



Now, using the constraint <p' 2 = —f, we obtain 



X((f>) = X e ™, 



(6.2) 



(6.3) 



(6.4) 



where Ao is an integration constant. For the above expression of \(<f>), junction conditions 
16.21 and 16.31 in fact, become independent of rj(t) (similar feature is reported in 



2. Junction conditions with Brans-Dicke field 



Similarly, in the Brans-Dicke case, the scalar field equation (Eq. 14. 15j) at the location of 
the brane (with the brane tension given by Eq. 16.11) gives us, 

H> 2 ] v dx 



a=ir b 



where 4> is, the mean value of the function (j)(t, a) at a = a b , defined as 



cr=a b 



Further, from the 00 or aa component of the Einstein equations, we get, 

3 [/I 

Using the relation 4l = nf, leads us to 

A(0) = Ao e~^\ 
for which the junction condition becomes, 

Ao 



[/'] = —q 7 ! e ~ 



(6.5) 



(6.6) 



(6.7) 



(6.8) 



(6.9) 



In this case, the junction conditions are time dependent. 

The junction conditions for the dilaton case can be obtained by using u = — 1 and a 
redefinition of the <ft field in the above analysis for the Brans-Dicke scalar. 
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B. Resolving bulk singularities 



It may be noted that all our solutions (except those discussed under the section on an 
exponential RS type warp factor) have bulk singularities, apart from the usual big-bang 
cosmological singularity in time. The latter is inevitable and expected while the former 
needs better understanding. 

The main question we need to address is - how do we manage to use a spacetime with a 
bulk singularity as a model for our bulk five dimensional geometry? There are ways to do 
this. 

One possible resolution of the singularities is done by putting branes at the locations of 
the singularities. For example, let us take, in the case of solutions with a bulk dilaton field, 
the two different solutions valid in two different regions (except the points where singularities 

llogk-^l+C, for a<a ^ 
— ~ log I cr - cr 2 | + c 2 , for a > a 
where <To is located somewhere between o\ and 2 . Then continuity at cr implies, 

^log |cr - oil + ci = -i log |cr - a 2 1 + c 2 , (6-11) 

Now this is a good bulk solution except at points 0\ and 2 . To resolve these singularities we 
can put two branes at o\ and a? and another brane at O and make their total contribution 
to bulk vacuum energy vanish 20] (as we do not have any bulk cosmological constant). This 
will be equivalent to putting delta-function source terms in the bulk Lagrangian. In effect, 
those singularities in bulk solutions actually provide us with places for the 3-branes, one of 
which can be chosen to be the Standard Model brane we live in. 



VII. CONCLUSIONS 



Finally, we list below, systematically, the conclusions obtained in this investigation. 

• We start out by writing down the energy condition inequalities using the Einstein 
tensor components (which have a nonzero Gq^). These inequalities are first checked 
for typical choices for the various metric functions a(t), r)(t) and f{&). We demonstrate 
that there exist viable models (following the criteria listed in the Introduction) which 
satisfy the energy conditions. 
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• We also find exact solutions with matter sources of various kinds. To begin with, we 
look at possible matter sources that may arise if we assume an exponential (in a) warp 
factor and some other typical constraints on matter stress energy. The solution space 
for [a(t),r)(t)] is analysed using a dynamical systems approach. Subsequently, we look 
at solutions with various types of scalar fields as sources - eg. bulk normal scalar, the 
Brans-Dicke scalar and the dilaton scalar. In our approach to obtaining solutions we 
exploit the known fact that with a conformally related metric and a redefined scalar one 
can convert BD theory into canonical Einstein Gravity coupled to a scalar. Analytic 
solutions are written down and the various possibilities that arise are outlined. In 
the case with a Brans-Dicke scalar and matter in the bulk, it is observed that, for 
traceless matter energy momentum exact and viable analytical solutions can be found 
with decaying or growing warp factors as well as decelerating/accelerating a(t). We 
note that in the case of the dilaton we can have decaying as well as growing warp 
factor solutions too. We are also able to find desirable solutions with growing a(t) and 
decaying r)(t) in several of our examples. 

• We check the energy conditions, the nature of the functions a(t) and r)(t) and figure 
out whether an obtained solution is desirable or not. The details are tabulated in one 
of the columns in the table commenting on the desirability of the solution. We do find 
several classes of desirable solutions, which satisfy our requirements and the energy 
conditions. 

• Finally, we outline how one can place a brane in the above bulk spacetimes by using 
the junction conditions . We also briefly discuss a way to resolve the bulk singularities 
using standard techniques. 

With three different functions appearing in the line element, it is always very difficult 
to find exact solutions. That we have found some is indeed encouraging. We hope to use 
these solutions in actual brane cosmological scenarios and arrive at relevant conclusions by 
making use of existing observational data on supernova and CMB anisotropies. We would 
also like to extend our results to other types of stress energy expressions. A disturbing 
aspect of our solutions is that they are singular in the bulk, though we do provide ways of 
resolving them. It remains to be seen whether we can find newer solutions (a) which are not 
necessarily of a power law type, (b) which have, in appropriate limits, both deceleration and 
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acceleration, (c) which have a decaying extra dimension stabilisable to a finite value and 
(d) where the bulk warp factor gives rise to a non-singular bulk metric, similar to RS. In 
essence, an appropriate combination of different bulk fields, with some dominating over the 
others in specific time intervals, will then be able to generate the expansion history (with 
proper decelerating and accelerating phase), of our universe on the brane. 
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